Magneto conductance oscillations in quasiballistic multimode nanowires 
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We calculate the conductance of quasi-one-dimensional nanowires with electronic states confined 
to a surface charge layer, in the presence of a uniform magnetic field. Two-terminal magnetocon- 
ductance (MC) between two leads deposited on the nanowire via tunnel barriers is dominated by 
density-of-states (DOS) singularities, when the leads are well apart. There is also a mesoscopic cor- 
rection due to a higher-order coherent tunneling between the leads for small lead separation. The 
corresponding MC structure depends on the interference between electron propagation via different 
channels connecting the leads, which in the simplest case, for the magnetic field along the wire axis, 
can be crudely characterized by relative winding numbers of paths enclosing the magnetic flux. In 
general, the MC oscillations are aperiodic, due to the Zeeman splitting, field misalignment with the 
wire axis, and a finite extent of electron distribution across the wire cross section, and are affected 
by spin-orbit coupling. The quantum-interference MC traces contain a wealth of information about 
the electronic structure of multichannel wires, which would be complimentary to the DOS mea- 
surements. We propose a four-terminal configuration to enhance the relative contribution of the 
higher-order tunneling processes and apply our results to realistic InAs nanowires carrying several 
quantum channels in the surface charge-accumulation layer. 



PACS numbers: 73.63.Nm,73.23.Ad,73.50. Jt 



I. INTRODUCTION 



Magnetoconductance (MC) of very small conductors, 
e.g., quantum dots, wires, and rings, is a topic of intense 
research since the early days of mesoscopic physics. In 
the simple case of phase-coherent conductors topologi- 
cally identical to the ring, one can generally distinguish 
between two classes of MC oscillations, differing in fre- 
quency by the factor of 2. First of all, the conductance os- 
cillates with a period of h/ e in magnetic flux (Aharonov- 
Bohm effect), since the magnetic field through the loop 
can be removed by a gauge transformation in integer mul- 
tiples of the flux quantum h/ e. These oscillations are 
sensitive to the disorder distribution, however, so that 
averaging over various disorder configurations in nom- 
inally identical conductors, or self-averaging in larger 
samples, reduces the strength of the oscillations. The 
quantum weak-localization (or antilocalization) correc- 
tion with a period of h/2e survives the averaging in dif- 
fusive samples.^ Depending on the concentration of the 
disorder in a given sample, either h/e- or /i/2e-type os- 
cillations can be more pronounced. A given realization 
of the disordered phase-coherent mesoscopic conductor, 
such as an open quantum dot or quasi-one-dimensional 
(ID) wire, also exhibits universal conductance fluctua- 



tions (UCF) with amplitude 



-/ft,, which vanish on 



sample averaging.-^. The weak-localization correction re- 
lated to h/2e MC oscillations in rings is usually mani- 
fested as a smooth background with the superimposed 
aperiodic UCF oscillations. 

All of the aforementioned phenomena can play a role in 



a quasi-lD wire with a finite cross section penetrated by 
a magnetic field. While a significant volume of the exist- 
ing literature is devoted to the MC studies of disordered 
wires (see, e.g., Refs. H and!3), the interference effects 
in magnetotransport of the quasiballistic nanowires did 
not attract the same attention (see, however, Refs. d for 
a discussion of weak-localization effects in narrow two- 
dimensional wires, which are ballistic across the wire 
width). Recent advances in nanofabrication techniques, 
however, make it possible to experimentally probe the 
MC of nearly ballistic semiconducting nanowiresr' rais- 
ing questions about the conductance oscillations related 
to the finite separation of the leads. See also Ref. 0, 
where ballistic carbon nanotubes were used to build a 
Fabry-Perot electron resonator, and Refs.H, reporting an 
observation and theory of the interference effects in tun- 
neling between parallel ballistic quantum wires through 
a finite-length barrier. 

In this paper, we are studying theoretically the mag- 
netotransport in single quasiballistic nanowires. Suppose 
the wire carries more than one quantum channel, so that 
the electron propagation between the injecting and de- 
tecting leads can interfere giving rise to conductance os- 
cillations in the magnetic field, which depend on the lead 
separation length. If the electron density is confined to 
a narrow surface layer and the Zeeman energy is dis- 
regarded, the MC in the parallel magnetic field can be 
decomposed into harmonics with frequencies in multiples 
of h/e approximately corresponding to relative winding 
numbers of various electron paths around the wire axis 
between the two leads. Generally, different transverse 
modes will effectively have different cross sections for ac- 
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cumulating the longitudinal magnetic flux resulting in 
aperiodic interference patterns. The aperiodicity is fur- 
ther enhanced by Zeeman energy and, furthermore, the 
spin degrees of freedom affect transport nontrivially due 
to the spin-orbit coupling. In addition, strong magnetic 
fields can have a significant effect on the one-dimensional 
band structure of a long wire, leading to van Hove mag- 
netofingerprints in the conductance. Varying the direc- 
tion of the magnetic field with respect to the wire axis 
gives an additional degree of freedom for probing elec- 
tronic structure. For example, a large perpendicular field 
eventually leads to a formation of Landau levels having 
a profound effect on the one-dimensional conductance. 

Among the principal practical questions are how the 
MC traces can be used to extract information about 
electronic states in the transverse direction (such as th 
e total number of modes and charge distribution, e.g., 
surface-confined or nearly uniform etc.) and about var- 
ious scattering mechanisms (involving in general elastic, 
phase-relaxation, and spin-orbit processes). While giv- 
ing comprehensive answers to these questions lies beyond 
the scope of this paper, our analysis can be used to con- 
sider them in some simple scenarios as well as to lay the 
grounds for more systematic studies. 

This work was motivated, in large part, by mag- 
netoconductance studies^ of cylindrical conducting 
nanowireSfS which provide a fascinating medium for 
mesoscopic physics^ i^°'^^ as well as show a great 
promise for applications as building blocks of nanowire 
networks. ^'^■^ After a general theoretical discussion in 
Sec. ini in Sec. IIIII we treat a specific case of trans- 
port in a surface charge-accumulation layer of cylin- 
drical InAs nanowires as those studied in Ref. _6, us- 
ing realistic parameters extracted from photoelectron 
spectroscopjii^ and magnetotransport measurements^^ 
on two-dimensional electron gases at InAs surfaces. Our 
approach may, however, also be useful for other types 
of conducting nanowires, such as state-of-the-art semi- 
conductor axial heterostructures with radially modulated 
composition and doping,^- Concluding remarks are given 
in Sec. IVl 



II. THEORY 

Consider a long ID wire connected through tunnel bar- 
riers to two narrow metallic source and drain leads sepa- 
rated by distance L, as shown in Fig. [T] In addition, the 
wire may be attached to two sinks (reservoirs) A and B, 
at the ends, whose purpose will become clear later. Disre- 
garding electron-electron interactions, the wire supports 
N transverse modes, the value of its two-terminal con- 
ductance being therefore bounded by Ne'^ /h, at a fixed 
magnetic field (supposing the sinks at the ends are elec- 
trically floating, for the moment). At a finite voltage bias, 
electrons are injected from the left lead 1 into the wire 
where they decompose into different transverse modes, 
are carried along the wire, and are detected by the right 



lead 2 or disappear in the reservoirs A and B at the ends 
of the wire. If the reservoirs are floating, an electron 
that enters the reservoir will eventually reenter the wire, 
but will have lost all coherence with the original injection 
process. The magnetic field applied either parallel or per- 
pendicular to the wire affects interference between differ- 
ent electron trajectories connecting the contacts, which 
accumulate the phase by enclosing a magnetic flux, re- 
sulting in conductance oscillations. In addition, the field 
can modify the number of quantum channels N at the 
Fermi level resulting in van Hove, i.e., density-of-states 
(DOS), singularities. 

We develop a simple formalism to capture interfer- 
ence between electrons injected into different transverse 
modes. Wave-function decomposition at the injecting 
lead is determined by the angular distribution around the 
wire circumference of the tunneling amplitude formed at 
the contact. (In a specific model discussed below, we 
take the tunneling strength to be uniform within an an- 
gle 29o determined by the lead deposition [see Fig.[Hb)].) 
For simplicity, we assume that the contacts are narrow 
along the wire in comparison to their separation L and 
to the electron wavelength, and we focus on MC oscilla- 
tions due to a finite length L. Similarly, the electrons are 
detected by the drain lead whose coupling to different 
transverse modes is determined by the deposition. We 
will focus on the MC oscillations due to the interference 
of phases accumulated during quasiballistic propagation 
between the leads. In particular, we will disregard any 
systematic contributions due to weak-localization (or an- 
tilocalization) physics, restricting our attention to oscil- 
lation features on scales determined by relative phases 
accumulated for multimode ballistic propagation along 
the length L. This physical picture is made mathemati- 
cally concrete in the following. 



A. General considerations 

A schematic of our model is shown in Fig. [T] A long 
wire of length and radius rg is contacted via tunnel 
barriers by two metallic leads 1 and 2 separated by dis- 
tance L and deposited on top of the wire [see Fig. [TJb)] . 
The total capacitance of the wire and/or the attached 
reservoirs is assumed large enough so that Coulomb- 
blockade effects can be neglected. However, L may be 
smaller than the coherence length, and we are interested 
in interference effects for electron propagation between 
contacts 1 and 2. 

Electron flow in the wire is carried by a surface charge 
layer. The disorder (due to impurities or wire-shape im- 
perfections) is assumed to be weak enough, so that the 
mean free path A > 27rro and it makes sense to define 
quasi-lD bands along the wire. If, furthermore, the lead 
separation L is not much longer than A, some electrons 
can propagate ballistically between the leads, giving rise 
to a magnetoresistance structure discussed in the follow- 
ing. 
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FIG. 1; Schematic of the model [view paraUel to the wire 
(a) and of its cross section (b)]: A long wire of length Lw is 
contacted via tunnel barrier by two leads 1 and 2 separated by 
distance L. In addition, the wire is 0hmically attached at the 
ends to two "sinks" A and B. Leads 1 and 2 can be voltage 
biased and A and B can be either grounded or left electrically 
floating. Sinks A and B are assumed to form good contacts 
with the wire and are used to either decohere (when floating) 
or sink (when grounded) incoming electrons, as explained in 
the text. 



We can distinguish between several regimes based on 
the relation between the contact separation L, mean- free 
path A, and the dephasing length l^. Let us assume that 
liji S> A, which can be achieved at low temperatures, and 
suppose the sinks A and B are floating. If L ^ 
A, the two-terminal conductance G12 is determined by 
summing three resistances in series (unaffected by A and 



^12 — + G2 



(1) 



where Gi is the tunneling conductance of the Zth lead and 
Gl oc 1/L is the diffusive conductance of the wire section 
between the leads. Gi can be experimentally obtained 
by measuring the two-terminal conductance between 1 
and A (while keeping the other contacts floating) and 
G2 — between 2 and B, for example, if the conductances 
of the sinks are large, Gi,G2 ^ Ga,Gb- The diffusive 
contributions from the wire can be accounted for by mak- 
ing additional two-terminal measurements between 1 and 
B and 2 and A, if the intercontact separations are known. 
(We are, however, interested in the clean limit, where the 
diffusive resistance of the wire is small.) If ^ L ^ A, 
Eq. ^ is still valid but the diffusive conductance Gl now 
acquires mesoscopic fluctuations that can be modulated 
by a nearby gate or the magnetic field. We will focus 
on the regime ^ X > L, where the electron propaga- 
tion between the contacts is phase coherent and nearly 
ballistic. 

To the lowest order in tunneling, the lead conductances 
Gl, I — 1,2, are determined by the density-of-states of the 
quasi- ID bands in the wire and the respective tunneling 
amplitudes Tn. For the ith ID mode in the wire, lead-^ 
tunneling conductance is Gu = \Tii\'^/vgi (with all the 
trivial prefactors including the lead DOS lumped into 



Tii), where 



(2) 



is the Fermi-point group velocity corresponding to the 
ith mode dispersion ei{k) (recall that the DOS is pro- 
portional to the inverse group velocity 1/vgi in ID). The 
tunneling conductance thus exhibits van Hove singulari- 
ties when the Fermi level crosses a band edge. If electrons 
are narrowly confined near the surface and are restricted 
to the lowest band in the radial direction, tunneling am- 
plitudes will be approximately the same for different ID 
bands and thus 



G, = 



(3) 



In this simple model, the two-terminal conductance G12 
given by Eq. ([1]) thus exhibits van Hove singularities 
whose magnetic fingerprints could be a useful tool for in- 
vestigating the electronic structure of the nanowires. Gl 
can also depend on the magnetic field, due to the weak- 
localization correction to semiclassical diffusion, but this 
will be disregarded assuming the wire is sufficiently clean. 

When the lead separation is decreased and becomes 
comparable to the mean-free path A, a mesoscopic cor- 
rection to the conductance SG12 , which is higher order in 
tunneling, may also contribute to the magnetoresistance 
traces for not too opaque tunnel barriers. 6G12 is deter- 
mined by a coherent tunneling between the leads, so it is 
fourth order in tunneling amplitude at the contacts and 
quadratic in the electron propagator in the wire. 



SG12 (X 



E 



(4) 



in terms of the retarded Green's function 



g,,{L,EF) = / die'^-*({vI',(L,t),*](0,0)}) (5) 
Jo 

for quasi- ID electron propagation between leads 1 
(source) and 2 (drain) at the Fermi energy Ep. (We 
set h — 1.) Mi^i^ is the tunneling matrix element into 
the ith channel in the wire from the lead-Z state labeled 
by /i at the Fermi energy. If the interband scattering in 
the wire is disregarded, Eq. @ becomes 



SG12 oc J2 



fiiy 



^Ml^^Mi,,,g,{L,EF) 



(6) 



At low bias and vanishing temperature, all quantities en- 
tering Eqs. (HI) and (O are evaluated at the Fermi energy 
since we are considering the elastic contribution to the 
tunneling. Although we will assume this in the following 
for simplicity, the generalization to the finite temperature 
and bias is straightforward for noninteracting electrons. 
(The temperature must, however, be at least larger than 



4 



the level spacing in the wire determined by L^, since 
we do not intend to analyze individual levels.) An addi- 
tional inelastic contribution to the higher-order tunneling 
(when one electron is injected into the wire but a differ- 
ent one is extracted, leaving behind an electron-hole pair 
excitation) may become important at finite temperature 
and/or bias,^^ but is disregarded here since it does not 
contribute to the discussed interference structure. For 
simplicity, we have assumed narrow leads, so that the 
propagator ([5]) is evaluated for a fixed distance L. The 
presence of sinks A and B at the wire ends decoheres 
electrons reflected from the ends and we can neglect their 
contribution to the propagator ([5]). 

If the incoherent contribution to tunneling determined 
by Eq. ([T]) is much larger than the correction (jj]) given 
by higher-order tunneling processes, the magnetoconduc- 
tance will be dominated by the van Hove singularities, 
according to Eq. ([3]) for the DOS-dependent contact con- 
ductances. It is however possible to suppress the inco- 
herent contribution to the signal by grounding sinks A 
and B at the ends of the wire and measuring the current 
in a grounded contact 1 or 2, while the other contact is 
voltage biased. The incoherent part of the conductance 
corresponding to the measured current will then be given 

by 

Q G1G2 ^ G1G2 

Gi + G2 -\- Ga + Gb Ga + Gb 

neglecting for simplicity the diffusive contribution of the 
wire and thus treating it as a disordered low-resistance 
node connecting four terminals. The approximation in 
the second equality is made assuming the leads 1 and 
2 to be much more resistive than the sinks: Gi,G2 ^ 
Ga,Gb- The coherent contribution jl]), on the other 
hand, is little affected by grounding the sinks since it is 
governed by the coherent electron propagation between 
the leads 1 and 2. Note that in the approximation of 
Eq. ([7]), G12 is proportional to the product G1G2 which 
is of the same order in the tunneling amplitude as the 
coherent correction SG12 given by Eq. Q. 

In the rest of the paper, we will mainly focus on calcu- 
lating the mesoscopic conductance ([6]) , which is sensitive 
to electron interference on traversing the distance L be- 
tween the leads. To that end, we will study the electronic 
structure in the long ID wire as a function of parallel or 
perpendicular magnetic field for certain simple models. 

In a clean ID wire, Qi{L,E) cx w^/, the inverse group 
velocity, so that the conductance ^ also appears sen- 
sitive to van Hove singularities. As we discuss in the 
following, however, disorder contributes to an exponen- 
tial decay of the propagator Gi{L, E) near the band edge, 
which suppresses DOS singularities in the conductance 
(|4]), making the interference effects more pronounced in 
the magnetoconductance traces. 

It is interesting to compare Eq. ^ to the result ob- 
tained for a different physical situation by Fisher and 
LeeJ^ These authors derive the Landauer-Biittiker ex- 



pression 

G^^X^I^^.f (8) 

for the conductance of a wire portion connecting two 
semi-infinite clean ID reservoirs, using Kubo formula. 
Here is the transmission-scattering amplitude between 
modes i and j. The transmission coefficient through a 
clean region ti oc VgiQi, for example, does not show any 
DOS singularities; in particular, |ip = 1 independently 
of the density for propagating modes, giving the cele- 
brated universal contribution to the conductance of /h 
per mode. The two-terminal conductance we are calcu- 
lating cannot, in principle, exceed the conductance ([8]) 
of the wire itself. The van Hove singularities obtained in 
the tunneling approximation are, in fact, artifacts of the 
perturbative expansion in the tunneling amplitude, and 
must saturate at a value below the wire conductance ([8]). 
Note incidentally that in the diffuse limit G^ corresponds 
to Gl introduced in Eq. ([1]). 

B. Cylindrical nanowires 

As a specific example, let us now consider a ballistic 
cylindrical wire with radius tq. Because of the transla- 
tional invariance and cylindrical symmetry of the wire, 
it is convenient to expand the one-electron wave function 
as 

0{r,e,z) = ^'^^Y.^V>nir), (9) 

71 ^ 

with a fixed momentum k along the wire axis z. (fnif) is 
the radial-position-dependent spinor with a given orbital 
angular momentum n along the wire. Assuming a uni- 
form magnetic field, we can express the translationally 
invariant Hamiltonian H according to 

Hip = e''^' Hn,n'^n' , (10) 

in terms of the 2x2 matrices -ff„.„' in spin space mixing 
orbital angular momenta n and n' . For a quasi- ID wire 
with the electron radial-confinement potential V{r), we 
obtain in a convenient gauge [Ag = Bzr/2 and A^ = 

2mr'^Hn^n' = [-r^df. + {krf + {n^4>,f - 1/4 + U{r) 
- 4'l{5n+2.ni + <5„-2,n' - 2(5„„') - mraRdrV{r) 

n,n' 

+ 2(j)xaxSn,n' - i'x (o-+(5„+2,n' + CT- (5„_2,n' )] 7 (H) 

where U{r) = 2mr'^V {r) is a dimensionless confining po- 
tential, (t>x,z{''') = e-Ba;^zr^/2 is the respective magnetic- 
field flux (in units of h/e) per area 7rr^ [thus note the 
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implied r dependence (j)x,z in Eq. (jlipj. — e is the elec- 
tron charge, rUf, is the free-electron mass, m is the effec- 
tive electron mass, g is the ^-factor, & = (iTa;, (Ty, CTz) is 
a vector of Pauli matrices, a± — ax ^ i^y, and B^^^) is 
the magnetic-field component parallel (perpendicular) to 
the wire. Notice that the perpendicular field Bx mixes 
angular-momentum states. The term proportional to 
Q!fl in Eq. (jlip is due to the Rashba spin-orbit cou- 
pling in cylindrically-symmetric confinement U{r) gov- 
erned by the material-specific coefficient an. Rashba in- 
teraction is believed to dominate spin-orbit coupling in 
InAs surface charge-accumulation layers due to a narrow 
gap and strong confinement fieldsi ^'^'^^i^^'^" The effective- 
mass approximation for the one-electron Hamiltonian 
()11|) breaks down near the wire surface, where the elec- 
tron confinement U{r) is sharp. Consequently, an appro- 
priate boundary condition for the wave function has to 
be applied at the surface (here, we assume ipn — at the 
boundary), while drU{r) defining the spin-orbit interac- 
tion is smooth within the electron-gas confinement. 

Depending on the confining potential V{r) and the 
wire radius, the radial electron-density distribution in the 
wire can be either peaked near the boundary or spread 
over the entire wire cross section. Figure [2] shows the 
zero-field electron-density profile within the wire, calcu- 
lated using parameters relevant to InAs nanowires (but 
for the purpose of demonstration setting ur — 0). The 
n = Q mode has some finite density p at the center of 
the wire [note that we are plotting p{r)r in the figure], 
while the finite-n modes do not contribute at the center 
and are pushed towards the boundary by the jlmr^ 
potential in the cylindrical coordinates [see Eq. pT| . 

The disorder-averaged Green's function for the ith 
mode is readily obtained after we calculate its dispersion 
ei(fc) from Eq. pT|) (assuming weak disorder), 

Ef - ei(k) + iT 

where the sum is over simple poles of [Ep — ei{k) + iF]^^ 
in the upper imaginary half plane [assuming no other sin- 
gularities under the analytic continuation of ei(fc) in the 
imaginary direction] and Vgi is the corresponding group 
velocity [see Eq. The disorder is accounted for in 

Eq. (|12p by an imaginary self-energy F = l/2r charac- 
terized by the scattering rate l/r. For weak disorder, 
we only take into account poles along the real axis with 
group velocity Vgi at ei{k) — Ep, the Fermi energy, and a 
corresponding imaginary contribution to the wave vector 
of 

Ski « iV/vgi . (13) 

In particular, we see that close to van Hove singularities 
where Vgi — > 0, the propagator (|12p is strongly damped 
due to disorder scattering. We will characterize the scat- 
tering strength by the mean free path A = vpT using 



15 




0.2 0.4 0.6 0.8 1 1.2 1.4 



FIG. 2: Electron-density distribution in the radial direction 
at zero magnetic field and vanishing spin-orbit interaction. 
The solid line is the total density p(r) multiplied by r (in 
arbitrary units) and the dashed lines are respective contribu- 
tions from the n — and n — ±1 orbital angular- momentum 
states. The Fermi energy Ef is chosen slightly below the 
n = ±2 band edge. The confining potential V{r) oc —r^^^ 
is flat at the origin and has the slope /3 near the hard-wall 
cutofi^ at r = ro- The arrow shows the radius containing half 
of the electrons. The specific choice of the parameters Ep, 
P, and ro is discussed in the beginning of Sec. IIIII re is a 
somewhat arbitrarily chosen "effective radius" of the electron 
distribution (see Sec. IIIII for more details). 



Fermi velocity vp of the two-dimensional (2D) electron 
gas corresponding to the surface charge layer. 

Since we base the rest of our discussion on the form 
(fT2)) of the disorder-averaged retarded Green's function, 
which determines the conductance 5Gi2 through Eq. ©, 
it is important to summarize our physical assumptions 
underlying these approximations. In Eq. ([S]), we only 
keep a coherent ballistic contribution that decays expo- 
nentially with intra- and intcrband scattering (thus ne- 
glecting vertex corrections). In particular, we disregard 
terms that are affected by multiple scattering between 
disordered regions. This is a major limitation if, e.g., 
there are strong-scattering regions at the contacts in- 
duced by the lead-deposition procedure. We have as- 
sumed that injected electrons that move to the left and 
those electrons that have passed under the right lead do 
not contribute to 6Gi2- Finally, cross-sectional shape 
and size variation of the wire could affect the interfer- 
ence pattern due to the magnetic field. For example, 
a strong radial variation can deteriorate the periodicity 
and strength of the magneto-oscillations in the parallel 
field if the electrons are confined at the surface. 

In summary, we study the conductance of translation- 
ally invariant nanowires, focusing on MC oscillations due 
to the ballistic paths between the leads. It is assumed 
that the paths that involve scattering from disorder result 
in rapidly oscillating contributions, which average out 
giving a nearly structureless background for the MC os- 
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cillations related to the baUistic interference (apart from 
the density-of-states singularities of the incoherent con- 
tribution Gi2 to the conductance discussed in the previ- 
ous subsection). If this is not the case, then our calcu- 
lations would apply to the mean conductance for a set 
of nominally identical samples. Electron-electron inter- 
actions arc not included, but are briefly commented on 
in Sec. HTdI 



C. Cylindrical-shell model 

As a further simplification, let us for a moment con- 
sider the shell model where electrons are confined in the 
radial direction to a narrow cylindrical well at radius 
re- Summing over the leads' degrees of freedom /i, 
in Eq. ^ and assuming the leads are metallic with high 
electron density, we arrive at the following expression for 
the coherent conductance: 

6Gi2 « J2g,g*{L,EF) f d9\Xs{e)\^0l{0)0j{e) 
X r de'\ue')\^\{e')^,{e'). (w) 



terms 



Here, (pi{9) is the angular component of the electron 
spinor wave function at the Fermi energy in the ith mode 
of the wire and \s{d){0) is the tunneling amplitude at the 
source (drain) lead, which is assumed to be independent 
of the electronic state in the leads but angle-dependent 
around the wire circumference. In the following, we set 
the tunneling amplitude to be Xs,d{0) = 1 for \9\ < 60 and 
zero otherwise, for both contacts. Finally, expressing the 
wave functions in terms of the angular-momentum eigen- 
states, (pi{0) = X^n'^inS*"^' obtain 



sin[(n — n')9Q] 



{n — n') 



If B,r 



(15) 

and the spin-orbit coupling is neglected, 
an — 0, the z component of the angular momen- 
tum commutes with the Hamiltonian (jlip . and there- 
fore Cin = \s)Sin, where |s) is the spin-s state along the 
z axis and Sin is the Kronecker delta. If we also ne- 
glect the Zeeman coupling, then by gauge invariance, 
the conductance of our shell-model wire must be a pe- 
riodic function of the flux (pz = so we may 
write 6Gi2i(t)z) = Z]m=-oo 5m exp(z27rTO0^). However, 
there may be considerable weight in the higher harmon- 
ics of the fundamental frequency, so that \gm\ may not 
be largest for small values of m. Also, the magnitude 
of SG12 may oscillate multiple times within the period 
— 1/2 < < 1/2, and it is natural to ask whether there 
may be a characteristic frequency of such oscillations. Al- 
though there is no simple answer to this question, it is 
helpful to look at the physical content of the interference 



sin[(r 



(16) 



between modes n and n' arising in this idealized situation 
(contributing to the conductance SG12 oc X^nn' G^"" ^). 
The last factor in Eq. (flB]) describes the suppression of 
the interference between modes with different orbital an- 
gular momenta due to the finite injection and detec- 
tion half-angle 9q. In particular, if 9q — n, all the 
cross terms vanish and G^"" ■* = |^nP<5nn', reminiscent of 
the resuli^ii for the Kubo conductance of a finite-length 
wire ([5]). Otherwise, we have to consider the oscilla- 
tory part of the cross terms proportional to ReC/„C/„': 

G(""') CX C0S[(A:„ - kn')L] exp[-rL(l/w„ + l/Vn')]/VnVn' , 

where kn is the Fermi wave vector and d„ group velocity 
for the nth mode along the wire. The corresponding rate 
of G*-"" •* phase accumulation in magnetic flux 0^ is given 
by 



UJ„ 



L d{kn - kn') 



27r 



d<Pz 



(17) 



For the shell model with weak scattering, fc^ — 
2mEF — (n + (j)z)'^/rl and, fixing Ep, (L/27r)(9fc„/902 — 
~{vg / Vn){L/2iTre), where mve = {n + (j)z)/re and mvn — 
kn- LOnn' is thus cqual to the relative winding number of 
the classical trajectories corresponding to modes n and 
n', on traversing the wire of length L around the cylin- 
der of radius re. Since w„„' is itself field dependent, it 
does not correspond to Fourier components of G*^"" ^ 
however. In addition, the latter is modulated by the mag- 
netic field via the prefactor exp[— rL(l/i;„-|-l/w„/)]/w„w„/ 
which has sharp features near the van Hove singularities 
in clean systems. Note that the fastest phase accumula- 
tion rates Unn' correspond to small velocities z)„ along the 
wire and are thus more sensitive to disorder scattering. 

We now consider the effect of a perpendicular field 
on the energy spectrum and tunneling DOS. We shall re- 
lax our assumption of a zero-thickness shell model, but 
we continue to assume that the confinement is sufficiently 
strong that at most one radial state is occupied for any 
given angular momentum n. Turning on the perpendicu- 
lar field Bx perturbatively gives the following correction 
to the energy eigenstate with angular momentum n (dis- 
regarding spin-orbit interaction): 



Ae„(fc) 



Am 



m 



n' — nil 



(18) 

which is readily obtained by second-order perturbation 
theory using Eq. ([TTjl [note that the lowest-order energy 
correction is quadratic in B^ since the term linear in 
Bx in the Hamiltonian (jlip mixes angular-momentum 
states that have different energies]. Here, (r^)„„ is the 
expectation value of the radial position squared in the 
nth state and (r)„„/ is the radial-position matrix el- 
ement between radial wave functions of states n and 
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n' . For example, applying this result to the shell model 

{{r^)nn = W)nn'? = ''e) gi^es (omitting Zeeman energy) 



20^ 



2mr'i 



2m 



(n + 0,)2-l/4 



where = eBxr^/2 is the flux that the perpendicular 
field Bx would produce in an area Trr^, in units of the 
flux quantum. We see that for a wire of circular cross 
section, a weak magnetic field in the perpendicular di- 
rection has less effect on the energies than the same field 
along the wire, since the former produces energy shifts 
proportional to 0^ while the latter produces shifts lin- 
ear in 02, for n ^ 0. The linear dependence on (j)z is a 
consequence of degeneracy between states of n and —n, 
which will in general be lifted by any deviations from cir- 
cular symmetry. The energy dependence on 0^ is then 
quadratic for 02 — > 0, but the curvature will be large if 
the zero-field splitting between n and —n is small, and 
the dependence on |02| becomes linear again when the 
shift due to 02 is larger than the zero-field splitting. 



D. Electron-electron interactions 

There can be several effects due to the electron-electron 
interactions affecting magnetoconductance oscillations in 
quasi-lD wires. In both disordered and clean wires, in- 
teractions lead to the inelastic-scattering rate growing 
with temperature. In the most naive form, this may be 
accounted for by adding a contribution to the scatter- 
ing rate F characterizing the single-electron propagator 
(fT2)) . A finite temperature can also lead to a qualitatively 
different effect, even in the absence of inelastic scatter- 
ing: The MC oscillation amplitude is suppressed as T"^/^ 
when ksT > i?Th, the Thouless energy, i.e., the inverse 
time for electron transfer between the contacts<^ 

For a wire that is connected to the outside world only 
through weak tunnel junctions, if the wire is short so 
that its Coulomb charging energy Ec — /2C is larger 
than the temperature, the usual Coulomb-blockade res- 
onant structure for the incoherent (sequential-tunneling) 
conductance G12 will set in as a function of source-drain 
and gate voltages. Away from the resonances, and at 
sufficiently low temperatures, higher-order quantum tun- 
neling of the electric charge (i.e., the cotunneling) will 
dominate the current, see, e.g., Ref. [H. For free elec- 
trons, the elastic component of the cotunneling signal 
reduces to our coherent conductance 5Gi2, Eq. The 
Coulomb-blockaded regime in short wires thus opens an 
interesting venue for studying quantum-interference ef- 
fects, which are higher-order in tunneling. 

At low temperatures and vanishing disorder, another 
interesting regime may be achieved characteristic of the 
nanowire geometry: that of the Luttinger-liquid state of 
the ID modesj^i In that case, the tunneling amplitude ac- 
quires a power-law suppression at each contact with low- 
ering temperature and bias,^^ determined by the anoma- 



lous compressibility of interacting electrons in ID. Some 
qualitative features of the ballistic magneto-oscillation 
patterns, however, may remain similar to the noninter- 
acting electrons. For example, the gauge invariance still 
dictates the h/ e periodicity in the parallel magnetic field 
for the narrow-confinement model, if we neglect the Zee- 
man energy. The spectral composition of the oscillations 
is, however, in general richer than discussed for nonin- 
teracting electrons in Sec. Ill CI In addition to the in- 
terference between one-electron trajectories between the 
contacts, there will also be contributions involving other 
electrons in the wire excited by the tunneling, which can 
similarly collect phase due to the applied magnetic field. 
Furthermore, the spin-charge separation characteristic of 
ID electron liquids may profoundly effect the interfer- 
ence at a finite voltage bias.* These topics are, however, 
beyond the scope of the present paper. 



III. NUMERICAL RESULTS 

As an application of the above formalism, we will here 
calculate the magnetoconductance in a model approxi- 
mating InAs nanowires. InAs surface (or interfaces with 
other materials) exhibits unusual properties, which are 
not manifested by most other semiconductors. InAs in- 
terfaces with metals typically do not form Schottky barri- 
ers, making it possible to achieve 0hmic contacts (which 
was utilized, e.g., in Ref. [l^ for using InAs nanowires 
as tunable supercurrent links between superconducting 
leads). Related to this is the property of InAs sur- 
faces to support two-dimensional quantum-well states 
whose properties can be tuned by the bulk doping levelJ^ 
Owing to this surface-charge layer, ID nanowires fab- 
ricated out of InAs can be expected to carry a finite 
number of ID quantum channels even when the bulk is 
insulatingi^i^ For example, taking the 2D electron density 
to be n2D = 10^^ cm.~'^r^^ we find that a 2D electron 
gas with radius Ve ~ 10 nm can support as many as 10 
quantum channels. 

In specific calculations, we will take the parameters'^ 
measured for a 2D inversion layer at the surface of a 
p-doped InAs bulk: Effective mass m — 0.027me, g- 
factor g — —15, Rashba coefficient au — 1 nm^, disorder- 
scattering mean-free path A — 100 nm, and triangular- 
well confinement with the slope drV — 15 meV/nm, 
which corresponds to the classical turning point of the 
ground state at the depth of about 10 nm. We will take 
the injection/detection half angle to be 9o = tt/5, lead 
separation L = 300 nm, and radius rp = 15 nm which, 
as explained below, corresponds to an effective radius of 
Te ~ 10 nm for electron distribution if using a quasi-2D 
shell model discussed in Sec. Ill CI (assuming that only 
the lowest quantum- well state is occupied). The char- 
acteristic spin-orbit splitting 2aiidrVkF is of the order 
of 10 meV. For simplicity, we will disregard the excited 
well states (which is reasonable for the chosen parame- 
ters), keeping in mind that the profile of the charge dis- 



8 




-6 -4 -2 2 4 6 



k[1/re] 

FIG. 3: Band structure of the ID wire in a large perpendicular 
magnetic field, (j>x ~ S (0z — 0), for the shell model with Ve — 
10 nm. Following a complicated nonmonotonic behavior of 
en{k) at small k and n, the energy grows nearly linearly with 
k at larger k, and eventually with the usual oc fc^ dispersion at 
k (X. These different regimes correspond to the dominance 
of the terms on the right-hand side of Eq. 1)11^ having the 
respective scaling with k. The dashed line shows the Fermi- 
level position. 

tribution can still be nontrivial on the scale of the wire 
radius. We fix the ID electron density, regardless of mag- 
netic field, to approximately correspond to a 2D electron 
gas with density n2D = 10^^ cm~^ residing at radius 
(more specifically, the density is chosen so that the zero- 
field lowest-band Fermi wave vector along the wire is the 
same as that of the 2D electron gas with density n2D)- 

A. Narrow- well confinement ("shell model") 

Let us first consider an idealized situation of a narrow- 
well confinement of the electron distribution at radius 
Te, which could be relevant for wider nanowires. As be- 
fore, we define the (dimensionless) magnetic flux to be 
(f'xiz) = ei?a;(z)rg/2. In the absence of the transverse 
field, (fix — 0, and neglecting the Zeeman splitting, gauge 
invariance dictates that the conductance is periodic in 
4>z with the period of unity. The periodicity is, how- 
ever, broken by a finite (7-factor and by an extended con- 
finement in the radial direction. According to the time- 
reversal symmetry, the linear conductance is symmetric 
under magnetic-field inversion, G{—B) = G{B), so that 
we need to consider positive fields only. 

The ID dispersion is given by Eq. (jl9p at low perpen- 
dicular fields, (t>x ^ and vanishing spin-orbit coupling, 
an = 0. In general, we can calculate the dispersion nu- 
merically by diagonalizing Hamiltonian pT|) (see Fig. [3] 
for — S). In order to find the conductance at a finite 
magnetic field, we first calculate the Fermi level keeping 
electron density fixed and then find the crossing points 
of en{k) = Ep with positive slopes. The Green's func- 
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FIG. 4: The main panel shows the interference MC traces in 
the parallel magnetic field for an idealized shell model with 
re = 10 nm, g = 0, and or = (solid line), a shell model 
including Zeeman splitting and spin-orbit coupling (dashed 
line), and a triangular- well model with outer wire radius 
ro = 15 nm (dotted line). All other parameters are listed 
in the text (in particular, L = 300 nm). The inset shows a 
calculation similar to the dotted line in the main panel but 
for the perpendicular magnetic field. The quantum-well states 
are treated in the WKB approximation (see Sec. IIII Bl for de- 
tails). 



tions for corresponding modes are then evaluated using 
Eq. (fT2l) and the conductance is finally obtained accord- 
ing to Eq. ()15|) , which requires knowing wave functions as 
well as dispersions [both given by the diagonalization of 
Eq. (fTTj) ]. Examples of the magnetoconductance traces 
are shown in Fig. [4] for parallel and perpendicular field 
orientations, taking into account both Zeeman splitting 
and spin-orbit coupling. In the case of the parallel field, 
the conductance is in general aperiodic in (j)^, although 
some reminiscence of the h/e periodicity is still present. 

In the Fig. m inset, there are considerable MC oscilla- 
tions at low perpendicular fields, before the conductance 
is suppressed at larger fields. The latter can be qualita- 
tively understood to be due to the onset of the Landau- 
level formation in the perpendicular magnetic field, which 
flattens the dispersions (see Fig. [3|), lowers group veloc- 
ity along the wire, and therefore, according to Eq. (fT^ . 
enhances disorder scattering on the propagation between 
the contacts. It is worthwhile noting in this regard one 
general feature near the van Hove singularities: Fixing 
the ID electron density, the Fermi level tends to be some- 
what attracted to the van Hove singularities before ID 
modes disappear or after they start getting populated, 
due to the diverging compressibility. The corresponding 
features can be seen in Fig. [5] in the Fermi- level position 
close to where it crosses transverse energy levels in the 
parallel magnetic field. 

As discussed in Sec. Ill Al see Eqs. ([!]) and the two- 
terminal source-drain magnetoconductance is expected 
to be dominated by the density-of-states singularities. It 
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FIG. 5: fc = energies for quasi- ID bands calculated by di- 
agonalizing Hamiltonian with the magnetic field along 
the wire axis for a shell model with = 10 nm. Dotted 
lines show spin-degenerate energies after setting g = and 
Or = 0. Solid lines are calculated taking into account the 
Zeeman and spin-orbit energies. Zero-field spin degeneracy 
is lifted by the latter while Zeeman splitting dominates at 
(j)z ^ 1- The dashed line is the Fermi energy correspond- 
ing to the solid-line energies, keeping the total ID electron 
density fixed. 
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FIG. 6: The total density of states as a function of the par- 
allel (black and dark-gray histograms) or perpendicular mag- 
netic field (light-gray histogram), which governs the incoher- 
ent two-terminal double-tunnel-barrier conductance ^ ac- 
cording to Eq. (Ul . The black histogram was calculated using 
the same parameters as the solid line in the main panel of 
Fig. m the dark-gray histogram — the dashed line, and the 
light-gray histogram — the dotted line. 



FIG. 7: Shell-model calculation of the interference magneto- 
conductance for 2D electron densities between 1 and 3 in 0.1 
increments (with lighter line shades corresponding to higher 
densities), in units of 10^^ cm~^. The dashed line is the aver- 
age. Here, g and an are set to zero, so that the 1 x 10^^ cm~^ 
trace is the same as the solid line in Fig. 

means, as a function of magnetic fields and B^, this 
would, in principle, allow us to obtain a wealth of infor- 
mation about the band structure and the electron dis- 
tribution in the radial direction [by using relations like 
Eq. (|18p obtained for low perpendicular fields B^]. A de- 
tailed procedure would require a further theoretical anal- 
ysis. 

Finally, before closing this subsection, we would like 
to note an interesting peculiarity of the magnetocon- 
ductance at low parallel magnetic fields. Varying the 
electron density, we find that there is more often a 
maximum of SGi2{4>z) sd, (j^z = 0; rather than a min- 
imum (there must be one of the two due to the time- 
reversal symmetry), over a wide range of parameters. 
This ballistic-interference maximum at low fields per- 
sists over density averaging, but it is different from the 
weak-antilocalization peak in the presence of a strong 
spin-orbit interaction in diffuse wires. See Fig. [7] for the 
shell-model calculation with zero spin-orbit coupling and 
Zeeman splitting, where the 2D electron density is varied 
between 1 and 3 in 0.1 increments, in units of 10^^ cm~^. 
We, however, expect the opposite statistical property, 
i.e., a local minimum at zero field for the incoherent con- 
tribution to the conductance G12. 



B. Triangular- well confinement 



is plotted in Fig. [HI By comparing to Fig. [H we see 
that the sharp DOS features are fully suppressed in the 
coherent-tunneling contribution to the conductance due 
to the significant disorder scattering. We note that if one 
was able to identify the positions of these van Hove sin- 
gularities in the two-terminal measurements or by other 



Consider now a triangular confinement potential 

V{r) = (3{ro - r) + ooe(r - tq) , (20) 

where (3 is the confinement steepness and is the Heav- 
iside step function defining the sharp cutoff at the outer 
radius tq. For the magnetic field along the wire axis, the 
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FIG. 8: Solid lines are the ground-state energies (relative to 
the n — energy) for transverse triangular-well states with 
fixed angular momenta n (labeled in the plot) as a function 
of the parallel magnetic field, setting g = and an = 0. The 
wire radius ro — 15 nm and other parameters are listed in the 
beginning of Sec. IIIII Dotted lines are calculated for the shell 
model with = 10 nm, which is a characteristic radius for the 
radial-distribution center of mass at small angular momenta 
n. See Fig. [2]for a similar confining potential V{r) (although 
not strictly triangular as here), (jiz = eBzr^/2. 

total effective potential in the radial direction depends on 
the angular-momentum state, the higher momenta (with 
respect to the flux — being pushed more outward, as 
can be seen from Eq. (fTTI) . We calculate the correspond- 
ing transverse ground-state energies for different angular 
momenta using the WKB approximation, in the absence 
of the spin-orbit interaction and perpendicular magnetic 
field. For simplicity, we use Eq. (|16p for suppression of 
the cross terms describing the interference between dif- 
ferent angular-momentum states, which was derived for 
the shell model. The spin-orbit coupling and, if present, 
the perpendicular magnetic field are also included sim- 
ilarly to the shell model, after we choose an appropri- 
ate effective radius r,, [by comparing triangular-well and 
shell- model energies (see Fig. HI)]. The calculated magne- 
toconductance for the parallel magnetic field is plotted 
in Fig. m by the dotted line. 

In Fig. [51 we plot the parallel-field dependence of the 
triangular-well ground-state energies at fixed angular mo- 
menta n, which were used in calculating the magneto- 
conductance. In order to compare WKB and shell-model 
calculations, the lowest energy is set to zero. For larger 
n (with respect to — 0z), the radial potential in Eq. (fTT|) 
pushes electrons outwards so that the electrons become 
confined to the narrow well near the surface at r = ro, as 
n ^ oo, decreasing the period of magneto-oscillations. 
This trend can be seen in Fig. [51 In particular, this 
means that it is rather crude to characterize different 
angular-momentum states by a single effective radius re 
and one would have to solve Eq. (|lip for a fully three- 
dimensional electron distribution (and beyond the WKB 



approximation, which is not expected to be very accu- 
rate for quantum-well ground-state energies in any case) 
if precise energy levels are desired. The zero-field radial 
electron distribution is shown in Fig. [21 for a somewhat 
more realistic confinement potential, which is fiat at the 
origin and normalized to have slope (3 at the boundary. 
The electron distribution is indeed rather broad within 
the wire. For a quantitative analysis, it would be nec- 
essary to determine the confining potential V{r) self- 
consistently with the electron distribution, as a function 
of the magnetic field. 

IV. DISCUSSION 

We theoretically studied the magnetoconductance 
(MC) of cylindrical quasi-one-dimensional nanowires car- 
rying several quantum channels confined near the surface. 
We focused on the MC oscillations arising due to the 
ballistic interference of electrons propagating via several 
ID modes along a finite-length wire section between a 
metallic source and drain leads. Van Hove singularities 
in the ID density of states are expected to dominate most 
prominent MC features in a two-terminal configuration, 
and we have therefore suggested to use additional reser- 
voirs at the ends of the wire to absorb injected electrons, 
which diffused away from the leads, in order to enhance 
the more interesting ballistic- interference features. 

In order to calculate mesoscopic MC interference pat- 
terns, we have studied the electronic structure of ID 
modes in the long wires for an arbitrary uniform mag- 
netic field, taking into account Rashba spin-orbit cou- 
pling due to the surface-confinement potential. In the 
case of the parallel magnetic field along the wire axis z, 
the MC is periodic in the magnetic fiux ipz with a period 
of h/e only in the case of a narrow surface confinement 
and vanishing Zeeman splitting. Relaxing these assump- 
tions, the conductance oscillations become aperiodic but 
some remnants of periodicity are still visible for realis- 
tic parameters, as seen in Fig. [4l A large perpendicular 
magnetic field leads to the Landau-level formation and 
suppression of the discussed mesoscopic magnetoconduc- 
tance. The amplitude of the conductance fiuctuations is 
reduced at a finite voltage bias and temperature or upon 
gate-voltage averaging modulating electron density. Cer- 
tain interference features, such as the local maximum in 
dGi2 at zero parallel magnetic field, however, can survive 
energy averaging (see Fig.[2l). 
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